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Higher-order vertices at zero external momenta for the scalar field theory, describ- 
ing the critical behaviour of the Ising model, are studied within the field-theoretical 
renormalization group (RG) approach in three dimensions. Dimensionless six-point 
ge and eight-point gs effective coupling constants are calculated in the three-loop ap- 
proximation. Their numerical values, universal at criticality, are estimated by means 
of the Fade and Pade-Borel summation of the RG expansions found and by putting 
the renormalized quartic coupling constant equal to its universal fixed-point value 
known from six-loop RG calculations. The values of gg obtained are compared with 
their analogs resulting from the e-expansion, Monte Carlo simulations, the Wegner- 
Houghton equations and the linked cluster expansion series. The field-theoretical 
estimates for g^ are shown to be in a good agreement with each other, differing 
considerably from the values given by other methods. 
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I. INTRODUCTION 



The critical thermodynamics of the three-dimensional Ising model is known to be de- 
scribed by Euclidean scalar field theory with the Hamiltonian 



where a bare mass squared uIq is proportional to T — Tj°\ Tj°^ being the phase transition 
temperature in the absence of the order parameter fiuctuations. Taking fiuctuations into 
account leads to renormalization of the mass ml — >■ m^, the field (pR, and the coupling 
constant A mg4, and also to the appearance of terms of the form 171^'"- g2n4''ji '^^^^ n > 2 in 
the effective action (free energy) of the system. In the critical region, where fiuctuations are 
so strong that they completely screen out the initial (bare) interaction, the behaviour of the 
system becomes universal and the dimensionless vertices g2n tend toward their asymptotic 
limits, i.e. they assume constant values which are also universal. 

In this paper, we calculate numerical values of Qq and gs at criticality using the field- 
theoretical renormalization group (RG) approach in three dimensions. Higher-order cou- 
pling constants mentioned related to corresponding vertices at zero external momenta will 
be found as power series in the renormalized dimensionless quartic coupling constant g^ up 
to three-loop order, and the RG series will be resummed by means of the Fade and Pade- 
Borel techniques. Then g^ in resummed RG expansions will be put equal to its universal 
value gl known from the canonical six-loop RG calculations resulting in numerical estimates 
for (yfg and g^. The numbers obtained for g^ will be compared with those resulting from the 
e~expansion for ge/gl as well as with the values given by Monte Carlo simulations, by the 
approximate solution of the Wegner-Houghton equations and by the analysis of the linked 
cluster expansions. 




(1) 
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II. RG EXPANSIONS FOR SEXTIC AND OCTIC EFFECTIVE INTERACTIONS 

Since higher-order bare coupUng constants are known to be irrelevant at criticahty in 
the RG sense, renormahzed perturbative expansions for qq, qs, etc. may be obtained from 
conventional Feynman-graph expansions of these quantities in terms of the only bare cou- 
pling constant - A. In its turn, A may be expressed perturbatively as a function of the 
renormahzed dimensionless quartic coupling constant ^'4. Substituting corresponding power 
series for A into original Feynman-graph expansions for qq and Qs, we can obtain the RG 
series for these higher-order effective coupling constants. 

Recently, one of the authors has found the sextic coupling constant Qq in the two-loop 
RG approximation [|l|] . Thus, what we are interesting in is the three-loop contribution to this 
quantity. As may be shown, there are 16 graphs contributing to in the three-loop order. 
In fact, their calculation is neither cumbersome nor lengthy since corresponding contribution 
may be written down as a sum of a few terms having a form of mass derivatives of some 
two-point and four-point graphs. So, within the three-loop approximation we get 



where go = A/4! and Z relates the dressed Green function G to the renormahzed one Gji in 
a conventional way: 



The renormahzed perturbative expansion for bare coupling constant g^ may be obtained 
using the normalizing condition 







go = mgi— 



(4) 



and the RG expansion for Z4 which has been calculated by many people: 




(5) 



Combining Eqs. 2, 4 and 5 we obtain 
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^6 = -9l(l - -94 + 1.389962952^1 



TT ^ TT 



(6) 



The RG expansion for the eight-point effective couphng constant gs may be found in an 
analogous way although it requires more job than in the case of ge- Two-loop and three- 
loop contributions to gs are given by 5 and 36 Feynman graphs respectively. Use of the trick 
mentioned above, however, considerably simplifies their calculation. The result is as follows 



With the RG expansions Eqs. 6, 8 in hand, we can get numerical estimates for universal 
critical values of ge and gg. 



Perturbative expansions in a field theory are known to be divergent, at best asymptotic. 
Moreover, the theory under consideration has no small parameter. That is why direct 
substitution of the fixed point value gl to Eqs. 6, 8 can not lead to satisfactory results. 
These expansions, however, contain important information which may be extracted provided 
some procedure making them convergent is applied. Here Pade and Pade-Borel methods 
will play roles of such procedures, i.e. we shall construct Pade approximants [L/M] for the 
functions given by the series Eqs. 6, 8 as well as for their Borel transforms which are related 
to the functions to be found ("sum of series") by the formula 





Expressing g^ in terms of g4 leads to the RG series for g^: 




(8) 



III. RESUMMATION AND NUMERICAL ESTIMATES 




(9) 



k=o ■ 



(10) 
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and then evaluate the integral Eq. 9 where series Eq. 10 possessing nonzero radii of conver- 
gence are replaced by corresponding Pade approximants. 

Starting from the three-loop expansions available, it is possible to construct Pade ap- 
proximants of the only reasonable type: [1/1]- On the other hand, one can try several 
different ways of resummation. Indeed, the Pade-Borel procedure may be applied not only 
to the series for qq and gs themselves but also to corresponding RG expansions for the ratios 
ge/gj and gs/g'l 

Let's construct first the Pade approximant for the six-point effective coupling constant. 
From Eq. 6 we readily obtain: 

9 ol + 0.500636^4 . ^ 

96 = -qi —■ (11) 

7r^M + 1.455566^4 

Substituting into this expression the most accurate numerical estimate gl = 0.988 for the 
fixed point value known from six-loop RG calculations in three dimensions [0,^ we find for 
the universal value of g^. 

gl = 1.694. (12) 

Apply then more sophisticated, Pade-Borel procedure to the series Eq. 6 which is ex- 
pected to lead to more accurate estimate for g^. After simple algebra we obtain 

o 7l + 0.052382o4t o , , 

c/6 = 0.477465^1 / — —t^e-'dt. (13) 

^ 1 + 0.291114^4^ 

Computation of this integral for g^ = 0.988 gives 

g; = 1.621. (14) 

The third way to get numerical estimate for g^ we use here is to construct the Pade-Borel 
approximant for ge/gf and to estimate the universal value of this ratio at criticality. The 
corresponding expression is as follows: 

gl 1 + 0.48519^4^ 



For the fixed point value of (74 tliis formula leads to: 

gl = 1.577. (16) 

To obtain numerical estimates for the same procedures are applied to the series Eq. 
8. They give 

g; = 0.68 (Fade), 

= 1.71 (Fade - Borel), (17) 

g; = 2.71 (Fade - Borel for gs/gl). 



IV. DISCUSSION 

Although the RG expansions found in Sec. 2 are, in fact, rather short, numerical estimates 
for g^ obtained from Eq. 6 by several resummation techniques are seen to be close to each 
other. That is why we believe that the numbers (12) and, in particular, (14) and (16) are 
also close to the exact value of g^. Moreover, taking into account of the three-loop RG 
contribution to g^ turns out to change corresponding estimate by less than 10 per cent 
(Fade-Borel resummation of the two-loop RG expansion have lead to g^ = 1.50 [|l|) what 
may be considered as an extra argument in favour of fair numerical accuracy of the results 
presented above. 

Let's compare our estimates for g^ with their counterparts obtained by different methods. 
The solution of the Wegner-Houghton equations within the local potential approximation 
presented by C.Bagnuls and C.Bervillier has yielded (ge/gl) = 3.59, gQ = 2.40 at the non- 
trivial fixed point . Determining by Monte Carlo simulations probability distributions for 
averaged magnetization of the 3D Ising model in an external magnetic field, M.M. Tsypin 
has found from his data that g^ = 2.05 The analysis of the linked cluster expansion 
series performed by T.Reisz has lead to g^ = 1.92 
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All these estimates are seen to be significantly larger than ours obtained by means of 
perturbative RG calculations for the field-theoretical model. In such a situation, it is in- 
teresting to evaluate using some alternative field-theoretical perturbative approach. The 
e-expansion for the ratio Qe/gl which results from the equation of state may be used for this 
purpose. Up to the order, it is as follows |]^: 

^ = 2e-^e2 + 1.2759el (18) 
gl 27 ^ ^ 

Resumming this expansion by Pade and Pade-Borel methods we find for e = 1 and g4 = 
0.988: 

g; = 1.687 (Pade), 

g; = 1.653 (Pade - Borel). (19) 

The Pade-Borel estimate for g^ thus obtained is remarkably close to those given by RG 
calculations in three dimensions (Eqs. 14 and 16). So, the field-theoretical RG approach in 
three and (4 — e) dimensions provides, within the three-loop approximation, the values of 
(jfg which agree well to each other but differ considerably from the Monte Carlo and some 
other estimates mentioned above. 

Let's discuss further the predictions concerning the eight-point effective coupling con- 
stant. The estimates following from the three-loop RG expansions in 3D are cleary seen 
to be strongly scattered being therefore numerically unreliable. It is not surprising since 
coefficients in the RG expansion for gs grow much faster than those in the RG series for 
gg making the former series less suitable for resummation than the latter. Correspondingly, 
much less reliable numerical estimates for gg at the critical point are obtained. On the other 
hand, all three numbers Eqs. 17 are i) positive, i.e. have a sign opposite to that of the 
lowest-order (one-loop) RG estimate g^ = —12.28, and ii) much smaller than the absolute 
value of this one-loop estimate. It means that higher-order RG contributions to gs being 
extremely important tend to strongly diminish its universal value at the critical point with 
respect to the number given by the lowest-order RG approximation. This conclusion seems 
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to be in accord with the fact that the recent analysis of relevant Monte Carlo data failed to 
reveal some appreciable (non-zero) value of 0. 

V. CONCLUSION 

In the paper, the RG series for coefficients before and in the expansion of the 
free energy of the three-dimensional Ising model in powers of the order parameter M have 
been calculated in the three-loop approximation. Numerical estimates for universal values 
of Qq and gs at the critical point have been found by Pade and Pade-Borel resummations 
of the series obtained and by putting the quartic coupling constant equal to its fixed point 
value 0.988. Pade-Borel procedure applied to the expansions for qq and ge/gf has given 
gQ = 1.621 and g^ = 1.577, respectively, while analogous treatment of the e-expansion for 
ge/gl has been shown to result in g^ = 1.653. Being in a good agreement to each other, these 
field-theoretical estimates are considerably smaller than those obtained by other methods. 
The RG expansion for g^ has turned out to be less suitable for resummation than that for 
gQ since it possesses rapidly growing coefficients. As a result, strongly scattered numerical 
estimates for g^ lying between 0.68 and 2.71 have been found. 
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